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Abstract 

We propose a mechanism that helps stabilize a superconducting state with broken time-reversal 
symmetry, which was predicted to realize in a d- wave superconducting film [A. B. Vorontsov, Phys. 
Rev. Lett. 102 , 177001 (2009)]. In this superconducting phase, the time-reversal symmetry break¬ 
ing is accompanied by spontaneous breaking of the translation symmetry along the film surface. 
We examine how the normal-superconducting phase boundary in the thickness-temperature phase 
diagram of the film is modified depending on the Fermi surface shape. In particular, the nonuni¬ 
form superconducting phase is found to substantially extend to a smaller thickness region in the 
phase diagram when the Fermi surface satisfies a nesting condition. We demonstrate this Fermi 
surface effect using a square-lattice tight-binding model. 
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Spontaneous time-reversal (TR) symmetry breaking in superconductors and superfluids 


& 


An 


was first established for the p + ip pairing state in the A phase of superfluid 3 He 
analogous chiral p -wave state has been discussed as a promising candidate for quasi-two- 
dimensional (Q2D) superconductor Sr 2 RuC >4 Q. The possibility of TR symmetry breaking 
is also discussed for heavy fermion superconductors PrOs 4 Sbi 2 js] and URu 2 Si 2 4] and for 


noncentrosymmetric superconductors LaNiC 2 [5| and Re 6 Zr [6J. 

Recent theoretical studies of the surface effects in unconventional pairing states have 
aroused renewed interest in TR symmetry-breaking states. Interestingly, the surface in su¬ 
perfluids and superconductors provides a mechanism responsible for spontaneous symmetry 
breaking. For example, when the B phase of superfluid 3 He is confined in the film geome¬ 
try, pair breaking at the film surface gives rise to spontaneous breaking of the translation 


symmetry along the surface |7]. Vorontsov found that a d -wave superconducting (SC) film 
with pair-breaking surfaces can exhibit superconductivity that breaks not only the trans¬ 
lation symmetry but also the TR symmetry js]. Hachiya et al. examined the stability of 
this TR symmetry-breaking phase against external magnetic fields jflj]. Very recently, a TR 
symmetry-breaking state accompanied by an unusual vortex pattern was predicted for a 


small d -wave SC grain 


10 ]. 


In this Rapid Communication, we address the effect of the Fermi surface shape on the 


phase transition of SC films. This study is motivated by previous theoretical studies 
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13] 


on the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) SC state that is stabilized in bulk materials 
under a strong applied magnetic field. The gap function in the FFLO state is reminiscent 


of those in nonuniform states proposed for the films of superfluid 3 He j?] and the d -wave 


superconductor 8]. One of the remarkable characteristics of the F 
upper critical field H c2 strongly depends on the Fermi surface shape 


'LO state is that its 


llMl3]. The question 


then naturally arises how the Fermi surface shape affects the phase transition in restricted 
geometries, where the instability is triggered not by an external field but by surface-pair 
breaking. We examine this problem in the context of the phase transition from the normal 


a. 


(N) state to the TR symmetry-breaking SC state predicted for a d -wave SC film 

The SC film considered here is the same as that in Ref. js] (Fig. [T]). Two parallel surfaces 
are located at y — ±D/2 and are assumed to be specular. The x-axis is taken along the 
surface, and the ^-direction is assumed to be uniform. We are interested in the phase 
boundary at which the N state becomes unstable against a Q2D SC state with the d -wave 
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FIG. 1: (Color online) A d -wave SC film with a distorted gap function. The gap function is strongly 
suppressed at the surface because of the destructive interference effect caused by quasiparticle 
scattering at the film surfaces. 


gap function 


A(R ,p) = 'y{p)A(R) 


( 1 ) 


with 

l{p) — V2sin2<£> p , (2) 

A (R) = A(y)e iqxX . (3) 

Here, R = (x, y) denotes the spatial coordinate, and <p p is the azimuthal angle specifying 
the direction of momentum p. The basis function y(p) is normalized as 

(7 2 (p)) = f Q y^7 2 b) = l- (4) 

The TR symmetry breaking in this state is due to the finite center-of-mass momentum q x 
of the d -wave Cooper pairs. 

Nagato and Nagai discussed the phase transition from the N state to the q x = 0 state in 
the him system of Fig. Q] with specular surfaces jwj]. At the N-SC phase boundary, A (y) 
was shown to take the form 


A(y) oc cos(q y y), q y = n/D. 


(5) 


The spatial variation of A (y) originates from the pair breaking caused by quasiparticle 
scattering at the him surfaces. As a consequence of the surface pair breaking, the him 
system has a critical thickness at which the N-SC phase transition occurs jl^|. Vorontsov 
:ound that the finite q x state has a smaller critical thickness than that of the q x — 0 state 
8|. This means that the N state instability hrst occurs for the q x ^ 0 state. 
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In previous theories fs[ Q, the Fermi surface is assumed to be isotropic (cylindrical). 
We generalize these theories to a him system with an anisotropic Fermi surface. To discuss 
the effect of the Fermi surface shape, we employ a square-lattice tight-binding model, which 


gives the dispersion 


3,Q, 


16| 


ip - —2t(cOSPz + cos Py) - p, 


( 6 ) 


where p is the chemical potential. In this model, cylindrical and square Fermi surfaces 
are obtained in the limits of p —> —4 1 and p —> 0, respectively. By controlling p, we can 
gradually change the Fermi surface shape from cylindrical to square (Fig. [5]). 

The N-SC phase boundary is determined from the gap equation 


T 


A(R) In — = tt 


x 


f(R,p,en ) - 


Po(<P P h(p) 

(po(^ P b 2 (p)> 

7 (p)A(il) 


(7) 


in which the pair amplitude f(R,p,e n ) obeys the linearized Eilenberger equation 17^191 


e n + -v{(p p ) • V R 


f{R,p,e n ) = sgn(e n ) 7 (p)A(i 2 ). 


( 8 ) 


Here, T c is the transition temperature in the bulk state, e n = 7rT(2n + 1) is the Matsubara 
frequency, and v(ip p ) is the Fermi velocity. Moreover, p 0 ((p p ) is the angle-dependent density 
of states at the Fermi level defined through the replacement 


y~l(' • •) j ^p('")^- (9) 

For the self-consistent solution A (R) oc cos {q y y)e iqxX , the linearized gap equation is 
reduced to 



where 


Vq = 2 V (^p) ■ Q = 2 M<Ap)<A + 


( 11 ) 


and A(<p p ) = p 0 (<p p ) 7 2 (p). Equation (TTOh determines the critical thickness as a function of 
(T, (p,,p). The p dependence comes from u(<p p ) and p 0 (<p p ). An optimum value of q x is 
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FIG. 2: Temperature dependence of the optimum q x for /i = —3.9 1, —0.3 1, and —O.Olf. 


determined such that the critical thickness is minimized at a given (T,/i), leading to a D(T) 
phase boundary for a given /i. 

In Fig. [21 we plot the optimum q x for /i = —3.9f, —0.3 1, and — O.Olf as a function of T/T c . 
The corresponding phase boundary lines are shown in Fig. [3j The vertical axis in Fig. [3] is 
7 t £ 0 /D = q y ^o with the coherence length dehned by 


£o = 


(^Op) + V y( ( Pp)) 

2nT, 


( 12 ) 


For fi = —3.9f, the Fermi surface shape is almost cylindrical, and the phase diagram [Fig. 
[3](a)] quantitatively agrees with that in Ref. js|. As the Fermi surface approaches the square 
shape, the phase boundary between the N state and the q x ^ 0 state (red solid lines) is 
remarkably enhanced, and simultaneously, the tricritical temperature T* is shifted higher. 
In the case of /i = —O.Olf, the phase boundary at T = 0 is shifted to 7 t^ 0 /D ~ 5.37. 

At T = 0, the critical thickness is determined by 


0 




(13) 


where 7 = 0.57721 • • • is Euler’s constant. In Fig. [4l we plot the critical value of tt£,o/D as a 
function of /i/t. As n/t approaches zero (the square-Fermi-surface limit), the critical inverse 
thickness increases steeply and diverges in the limit ji/t —> 0 . 

The reason for phase boundary enhancement due to the anisotropic Fermi surface can be 
understood as follows. In the him system under consideration, the gap function is suppressed 
by surface scattering. In Eq. (ITU]) , the surface pair-breaking effect is described through r] q 
with q x = 0 , i.e., 


Vq y = v y ^Pp)%l c 2 = Vy((p p )ir/2D. 


(14) 
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FIG. 3: (Color online) Phase diagram for (a) fj, = —3.9f, (b) // = —0.3 1, and (c) fx = —O.Olt. The 
red solid (blue dashed) curve is the phase boundary between the N state and the q x / 0 (q x = 0) 
state. The arrows indicate the tricritical temperature T* « (a) 0.43T C , (b) 0.57 T c , and (c) 0.59T C . 
In the inset, we depict the Fermi surface corresponding to the /r value. 

The finite modulation q x adds the term r) qx = v x (tp p )q x / 2 to rj q . This additional term can 
cause r] q to equal zero on some portions of the Fermi surface. The condensation energy lost 
by the surface pair breaking can thus be compensated by inducing modulation along the 
surface. This accounts for the stabilization of the finite q x state in the film. Note that the 
condition i] q = 0 on the Fermi surface is equivalent to the nesting condition £ p — £= 0 
of the Fermi surface. Thus, good Fermi surface nesting helps stabilize the finite q x state. In 
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FIG. 4: Critical value of tt^q/D at T = 0 as a function of n/t. 
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FIG. 5: (Color online) Fermi surface nesting for fi = —3.9 1, —2 1, —0.3 1, and —O.Olt. The red dotted 
curves show the Fermi surface shifted by q = (q x , q y ). 

Fig.0 we show the Fermi surfaces for several // values. The nesting condition can be satisfied 
in a finite area for /i = —O.Olt, whereas it can only be satisfied on lines for /j = —3.9 1 as 
in the case of the isotropic (cylindrical) Fermi surface. As a result, the nearly square Fermi 
surface (/i = —O.Olt) causes substantial enhancement of the phase boundary between the N 
state and the q x ^ 0 state. 
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A similar Fermi surface effect has been discussed in the context of FFLO instability 


111 Il3| . In that case, the surface pair-breaking term r) q in i] q is replaced by the Zeeman 


coupling to an external magnetic field, and the critical inverse thickness corresponds to the 


upper critical field H c2 . As shown with the FFLO problem for Q1D systems 


20H22|, the 


H c2 (T) curve has a positive curvature in the case of perfect nesting. The corresponding 
behavior is found in Fig. |2](c) (/i = —O.Olt) in the upturn of the red solid line below T*. 
When /i = —0.3t, similar upturn behavior is found; however, in this case, the curvature 
becomes negative at low temperatures. The low-temperature difference occurs because the 
nesting condition for /i = —0.3t is not “touching on surfaces” but “crossing on lines.” At 
high temperatures, the thermal energy k^T makes the difference between the two conditions 
indistinguishable [13|, and the phase boundary line exhibits an upturn similar to the case of 
fi = —O.Olt. 

In conclusion, we have discussed the effect of the Fermi surface shape on the N-SC phase 
transition in unconventional SC films with pair-breaking surfaces. We have demonstrated 
the Fermi surface effect using the phase transition to the TR symmetry breaking state 
predicted for a d -wave SC him 8] as an example, and we showed that the critical thickness 
is substantially reduced near half hlling, where the Fermi surface is almost square shaped. 
This result can be interpreted as a consequence of Fermi surface nesting. An analogous 


Fermi surface effect has been discussed for the FFLO stability problem 


11 


13]. In this case, 


Fermi surface nesting causes an enhancement of the upper critical held H c2 . In both cases, 
the spatial modulation of the gap function is induced to avoid strong pair breaking. The 
pair breaking is caused by surface scattering in the him case or by an external magnetic 
held in the FFLO case. In general, Fermi surface nesting allows a much greater reduction 
in the pair-breaking effect. The stabilization scenario for nonuniform superconductivity by 
the Fermi surface nesting can be applied quite generally to superconductors with strong pair 
breaking. 
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